High-Resolution Magnetometry with a Spinor Bose-Einstein Condensate 
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We demonstrate a precision magnetic microscope based on direct imaging of the Larmor preces- 
sion of a *^Rb spinor Bose-Einstein condensate. This magnetometer attains a field sensitivity of 
8.3 pT/Hz^/^ over a measurement area of 120 /^m^, an improvement over the low-frequency field 
sensitivity of modern SQUID magnetometers. The corresponding atom shot-noise limited sensitivity 
is estimated to be 0.15 pT/Hz^/^ for unity duty cycle measurement. The achieved phase sensitiv- 
ity is close to the atom shot-noise limit suggesting possibilities of spatially resolved spin-squeezed 
magnetometry. This magnetometer marks a significant application of degenerate atomic gases to 
metrology. 



Precision magnetometers that map magnetic fields 
with high spatial resolution have been applied to studies 
of condensed matter systems biomagnetic imaging 
and tests of fundamental symmetries Many of 
these applications require the measurement of magnetic 
fields at low (< 10 Hz) frequencies. Current technologies 
capable of micron-scale magnetic microscopy include su- 
perconducting quantum interference devices (SQUIDs), 
scanning Hall probe microscopes, magnetic force micro- 
scopes and magneto-optical imaging techniques Of 
these, SQUIDs offer the highest sensitivity, demonstrated 
at 30 pT/Hz^/^ over a measurement area of around 100 
/im^ £| . The low-frequency sensitivity of these devices is 
limited by (1//) flicker noise of unknown origins 6]. 

Magnetic fields may also be sensed by detecting the 
Larmor precession of spin polarized atomic gases. To 
date, atomic magnetometers have achieved field sensitiv- 
ities of 0.5 fT/Hz^/^ over measurement volumes of 0.3 
cm'^ 0. However, attaining high spatial resolution with 
a hot-vapor medium is precluded by rapid thermal dif- 
fusion of the atoms, restricting the minimum resolved 
length scale of these magnetometers to around 1 mm. 

Trapped ultracold gases present an attractive medium 
for a variety of precision measurements due to their 
negligible Doppler broadening and long coherence times 
Il3|- Spinor Bose gases, comprised of atoms with 
non-zero spin, the orientation of which is free to vary, are 
particularly well suited to magnetic microscopy. In con- 
trast with hot- vapor atomic magnetometers, the suppres- 
sion of thermal diffusion in a gas through Bose-Einstein 
condensation enables precise measurements at high spa- 
tial resolution. Also, density-dependent mean field shifts, 
which deleteriously affect other types of precision mea- 
surements using dense ultracold gases, do not affect Lar- 
mor precession due to the rotational invariance of inter- 
particle interactions in a spinor gas 
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Here, we perform precise magnetic microscopy with 
high two-dimensional spatial resolution using a ^'^Rb 
F = 1 spinor Bose-Einstein condensate (BEC). In our 
magnetometer, longitudinally spin-polarized spinor con- 
densates are prepared in an optical trap. Larmor pre- 



cession is induced by tipping the magnetization perpen- 
dicular to a bias field imposed along the axis of the con- 
densate. The spins in each region of the condensate then 
precess at a rate that is proportional to the local mag- 
netic field. After a variable integration time, the conden- 
sate is probed using magnetization-sensitive imaging to 
extract the local Larmor phase. The difference in this 
phase between various regions of the condensate reveals 
the spatial variations of the magnetic field. 

The determination of the accrued Larmor phase of a 
coherent spin state, such as the transversely magnetized 
condensate, is subject to an uncertainty in the initial 
phase of S(f>a — l/ViV due to projection noise of measur- 
ing N atoms. This noise limits the field sensitivity over a 



measurement area A to 6B 



where r is 
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the Zeeman coherence time and h the local column den- 
sity of the gas. We assume the measurement is repeated 
over a total measurement time T at a duty cycle D. The 
scaling of field sensitivity with the measurement 
area A for the atomic magnetometer may be compared 
with the area scaling for SQUID magnetometers. This 
scaling ranges between A'^/^, for a fixed SQUID sensor 
coupled optimally to a variable pickup loop, and A~^/^, 
for direct sensing with a SQUID optimized to operate at 
the quantum limit for the noise energy [l3| . For either 
scaling, the atomic magnetometer outperforms SQUID 
magnetometers at small measurement areas (Fig. [1]). 

Optical detection of Larmor precession is limited also 
by photon shot noise. In this work, the Larmor precession 
phase is measured by repeated phase-contrast imaging of 
the condensate using circular polarized light 0] . For our 
probe detuning of (5 = 27r x 500 MHz below the F = 
1 F' — 2 (Dl) transition of ^^Rb, the phase contrast 
signal can be written as s ~ 1 -I- 2hcro{'y /2S)[co + ci{Fy)] 
where co = 3A^/27r is the resonant cross-section, 7 is the 
natural linewidth and Fy is the projection of the local 
atomic spin on the imaging axis y, which is made perpen- 
dicular to the field axis. The detuning-dependent con- 
stants Co — 0.118 and ci = 0.274 describe the isotropic 
polarizability and optical activity, respectively. We ne- 
glect the effects of linear birefringence (ex {Fy)). An es- 
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FIG. 1: Field sensitivity for repeated measurements using 
the spinor BEG magnetometer. Gurves marked ASN (PSN) 
represent atom (photon) shot-noise limited sensitivities, as- 
suming r = 250 ms, D = 1, and the atomic column density 
and probe light levels for our experiment (see text). Diffu- 
sion of magnetization limits the sensitivity for a given length 
scale by imposing a limit on r (short dashed line, assuming 
D = 1). The gray line indicates the measured spatial root Al- 
lan variance; the sensitivity demonstrated in measurements 
of an applied localized magnetic field, assuming duty cycles 
oi D = 0.003 (•) or D — 1 (o), is also shown. Results are 
compared both to the ideal sensitivity of a quantum-limited 
SQUID magnetometer (dot-dashed line) and to demonstrated 
low- frequency sensitivities [1,[3 ('^)- 

timate of the Larmor precession phase is obtained by 
tracking the sinusoidal oscillation of the phase contrast 
signal across the sequence of phase contrast images. The 
photon shot-noisc limited sensitivity of this estimate is 
then (50^ w ^_ v__^_^_^_^ gjyij^g a field sensi- 
tivity limit of SB — — — ££2=_ Here, ri is the detection 
quantum efficiency and Np is the total photon fluence, in- 
tegrated across the multi-pulse imaging sequence, vifithin 
the region of interest. 

For our demonstration, spin-polarized *''Rb conden- 
sates of up to 1.4 X 10^ atoms were confined in a single- 
beam optical dipole trap characterized by trap frequen- 
cies {uj^,ujy,uj^) = 27r(165,440,4.4) s'^ 0. The tight 
confinement along the imaging axis (condensate radius 
Vy = 2.0 fim) ensured that the condensate is effectively 
two-dimensional with respect to spin dynamics. Larmor 
precession of the condensate was induced in the presence 
a bias field of 165(7) mG aligned along the long axis of 
the condensate. A measurement integration time of 250 
ms was chosen; at longer times, measurements were ham- 
pered by uncontrolled motion of the condensate along the 
weakly confining dimension (see below). 

We operated our ultracold-atom magnetometer under 
two testing conditions. In one, the magnetometer was 
used to sense the long length-scale inhomogeneous back- 
ground magnetic field in our apparatus. While these 
background fields were partially cancelled by suitable 
electromagnets, the remaining background typically var- 



ied from shot to shot. To account for this fluctuating 
background, the field profile along the long axis of the 
condensate was determined by a third-order polynomial 
fit to the magnetometer measurements [isj. The residu- 
als from this fit were then analyzed to characterize noise 
limits to our magnetometer. 

In the second testing condition, we used the magne- 
tometer to measure a deliberately applied, localized mag- 
netic field. This field was simulated using a circularly 
polarized laser beam at a wavelength of 790 nm. The 
choice of wavelength and polarization ensured that this 
beam imposed a local optically-induced Zeeman shift 
on the trapped atoms (Fig. [2^). The beam was aligned 
at an angle 9 ~ 60° to the direction of the bias field, 
incident and focussed in the plane perpendicular to the 
imaging axis. The magnetic background for each run of 
the magnetometer was again determined by third-order 
polynomial fits, but using measurements from regions far 
from the localized field. The Zeeman shift due to the lo- 
calized field was extracted from the residual of this fit. 

Measurements of this localized field were affected by 
small center-of-mass oscillations of the condensate along 
its long axis. An oscillation with amplitude Sz blurs the 
magnetic landscape and washes out features comparable 
to or smaller than 5z. Unable to eliminate such oscil- 
lations completely, we monitored the condensate motion 
for each run of the magnetometer by a sequence of four 
images spaced at a quarter period of the axial trap fre- 
quency, taken prior to the Larmor imaging sequence. We 
discarded measurements for which an excursion compa- 
rable to the width of the localized field was indicated. 

Two-dimensional maps of the magnetic field were ob- 
tained from a pixel-by-pixel analysis of the Larmor pre- 
cession phase at each coordinate within the profile of the 
condensate (Fig. [2)3). The frame-to- frame variation of 
this signal showed the characteristic oscillations due to 
Larmor precession as well as an overall decay of the con- 
densate number due to off-resonant scattering of probe 
light. This decay was taken into account in obtaining an 
unbiased estimate of the local Larmor phase. Our 2D 
approach was found to be susceptible to imaging aberra- 
tions, primarily in the narrow (x) dimension of the gas. 

Alternately, more robust measurements were obtained 
by integrating the field measurement over the x direc- 
tion to reduce measurements to a single resolved direction 
along the long axis (z) of the gas. For this, the aberrated 
signal profile in the x direction was determined at each z 
coordinate in the images from averages over the multiple 
frames. The phase-contrast signal height in each image 
frame and at each z coordinate was then determined. ID 
phase profiles were obtained as before (Fig. [2t) . 

The demonstrated sensitivity of our magnetometer is 
shown in Fig. [TJ The spatial root Allan variance from 
the ID data was determined for each of 15 runs of the 
magnetometer under the first testing conditions (mag- 
netic background only) and then averaged. Here, the 
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FIG. 2: (a) Sequence of phase contrast images of a condensate 
under the influence of a local optically-induced Zeeman shift, 
(b) The resulting 2-D map of the magnetic fleld obtained by a 
pixel- by-pixel estimation of the Larmor precession phase, (c) 
1-D phase profile showing the applied field inhomogeneity. 
The peak strength of the field as determined by this single- 
shot measurement is 166.2 ± 1.2 pT. 



measurement area is determined by accounting for the 
effective 5.3 /xm length over which the aberrated signals 
are averaged in the x direction. This noise level matches 
closely with photon shot-noise estimates and is ^--^3 times 
that due to atomic shot-noise given the number of atoms 
in the corresponding areas. Excess noise for areas larger 
than about 20 /itm^ was found to correlate with the lo- 
cal intensity of the probe light, an effect we attribute to 
probe-light induced shifts of the Larmor precession fre- 
quency during imaging. This noise can be reduced by us- 
ing a linearly polarized probe with a more homogenous 
intensity profile and by carefully aligning the magnetic 
bias field to be perpendicular to the imaging axis. 

Results from measurements under the second testing 
condition (background plus localized field) are shown in 
Fig- [3 Here, the strength of the applied field (peak value 
of Gaussian fits) was measured repeatedly at several pow- 
ers of the field-inducing laser beam. From all these mea- 
surements, a calibration between the laser power and the 
localized field strength was obtained. From the residual 
scatter in the measurements, we determine the (rms) sen- 
sitivity of our Larmor precession phase measurements as 
1.0 X 10~^ rad, corresponding to a single-shot field sensi- 
tivity of 0.9 pT over the 120 ^m^ area under the Gaussian 
field profile. This sensitivity was demonstrated for field 
strengths up to 60 pT. A marginally larger variance for 
higher fields points to the existence of small systematic 
effects, e.g. residual motion of the condensate or varia- 
tions in the localized field strength. Averaging over the 
entire range of measurements shown in Fig. [3] yields a 
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FIG. 3: Single-shot measurements of the localized magnetic 
field imposed by a laser beam focused to cr = 24 fim rms 
width. Residuals from a linear fit, calibrating the field magni- 
tude vs laser power, are shown. Error bars indicate standard 
deviations for 10 measurements per setting. The measure- 
ment time was 250 ms. 



phase sensitivity of 1.2 x 10 ^ rad. 

Under repeated operation, our magnetometer, with a 
low duty cycle of just Z) = 3 x 10^^, attains a field sensi- 
tivity of 8.3 pT/Hz^/^, an improvement over that demon- 
strated for low-frequency (< 10 Hz) field measurements 
with modern SQUID magnetometers 0, [l3|. Plausible 
extensions of current cold-atom experimental methods 
should enable duty cycles of order unity. At full duty cy- 
cle, our demonstrated single-shot sensitivity would cor- 
respond to a field sensitivity of 0.5 pT/Hz^/^. 

As discussed above, in the photon shot-noise limit, the 
sensitivity of an atomic magnetometer increases with in- 
creasing probe fluence. While calculations based on lin- 
ear Raman scattering rates indicated that reliable phase 
measurements could be obtained even at a fluence of 3400 
photons//xm^ per frame, it was found that the light in- 
duced losses of the condensate far exceeded those pre- 
dicted by the calculations. The discrepancy was at- 
tributed to superradiant Raman scattering of atoms into 
the F — 2 hyperfine states, in which atoms would no 
longer be observed by our probe. To counter this prob- 
lem, the superradiant gain was reduced by lowering the 
probe intensity. Each frame in the imaging sequence 
was obtained by integrating the light from four pulses 
of light, each of duration 2.2 /is and spaced by the Lar- 
mor period of ~10 /iS. We also increased the motional 
decoherence of the F = 2 atoms produced during super- 
radiance by their preferential scattering of incident light 
resonant with the F = 2 ^ F' = 3 (D2) transition. To- 
gether, these strategies enabled a probe fluence of 750 
photons//im^ per frame. 

Our magnetometry medium, though Bose condensed, 
is still a gas in which atoms are free to move. Thus, in 
determining the phase shift accrued due to a local mag- 
netic field, one must consider atomic motion due to both 
quantum-mechanical and classical effects. For instance, 
imposing a weak inhomogeneous field of characteristic 
length a leads to quantum diffusive motion of the fluid. 
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For times t > tq = ma'^ /h, with m the atomic mass, the 
motion of the spinor gas will greatly reduce the phase 
accrued due to Larmor precession. This evolution can be 
considered to be the quantum limit of thermal diffusion 
observed in NMR studies 17|. In our experiment, the 
condition r = 250 ms > tq is reached at length scales 
below 10 ^m (Fig. [1]). For integration times t tq, 
effects of quantum diffusion require that phase measure- 
ments (f> be corrected by an amount (5(/>/0 ~ {'''/tq)'^- 

An inhomogeneous magnetic field also exerts forces on 
magnetic dipoles. In the extreme case of static inho- 
mogenous fields with a Zeeman energy comparable to 
the chemical potential (~ 10^ times larger than those 
studied in this work), these classical forces can modify 
the density distribution of the condensate [l^. In our 
case, these forces result in small corrections and limits 
to the accrued Larmor phase. The field strength B and 
a characteristic length cr for its variation define a classi- 
cal time scale tc = \/ ma'^ / hbB , the time taken by an 
atom to move a when accelerated by this field. For an 
integration time r, this classical motion imposes a limit 
on the maximum detectable phase shift (when t = tc) 
of (jim — tq/t. It should be noted that neither the dif- 
fusion of an imprinted phase nor the limitation on the 
dynamic range are fundamental; both can be eliminated 
by constraining atomic motion, e.g. by imposing an op- 
tical lattice potential. 

To observe the dilution of magnetization due to atomic 
motion, we imposed a light-induced magnetic field local- 
ized to a beam waist of 5.4 /im onto the transversely- 
magnetized spinor condensate. Following a 5 ms ex- 
posure to the field-inducing laser beam, the condensate 
magnetization was allowed to evolve freely for variable 
time before being probed. During this evolution, the 
imprinted Larmor phase diminished in peak height and 
grew in extent, matching well with calculations based on 
a non- interacting spinor gas in a localized field (Fig. 2]) . 

In conclusion, we have demonstrated a spinor-BEC 
magnetometer, a powerful application of ultracold atoms 
to precision measurement of scientific and technological 
significance. Inasmuch as the Larmor precession phase 
represents the phase relations among BECs in several 
Zeeman states, this magnetometer can be regarded as a 
condensate interferometer with high temporal and spatial 
resolution. The single-shot phase sensitivity and shot- 
to-shot variations of 10 mrad achieved here represent an 
order of magnitude improvement over the performance of 
current EEC interferometers. 

The demonstrated phase sensitivity is close to the atom 
shot-noise limit. This augurs spin-squeezed magnetom- 
etry via continuous quantum non-demolition measure- 
ments of the condensate [19] and novel spatially and tem- 
porally resolved studies of spin squeezed ensembles. Fi- 
nally, the high spatial resolution and sensitivity of the 
spinor magnetometer render it a powerful detector of 
quantum spins in optical lattices and may facilitate time- 
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FIG. 4: Quantum evolution of an imprinted phase, (a) Cross 
sections of the phase imprint after a free evolution time t = 
0,40 and 90 ms. Traces are offset for clarity, (b) The 1/e 
width of the imprinted phase. The dashed line corresponds 
to numerical simulations based on a non-interacting spinor 
gas. (c) The peak phase vs the free evolution time. The 
dashed line corresponds to the expected scaling imposed by 
normalization of magnetization. 



resolved, nondestructive studies of frustration, dipolar in- 
teractions and disorder in such systems. 

We thank D. Budker and J. Clarke for helpful discus- 
sions. This work was supported by the NSF and the 
David and Lucile Packard Foundation. S. R. L. acknowl- 
edges support from the NSERC. 



[1] 

[2] 

[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 
[12] 

[13] 

[14] 

[15] 



[16] 
[17] 



These authors contributed equally to this work. 

C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 

(2000). 

K. Kobayashi and Y. Uchikawa, IEEE Trans. Mag. 39, 
3378 (2003). 

L. R. Hunter, Science 252, 73 (1991). 
S. J. Bending, Adv. Phys. 48, 449 (1999). 
J. R. Kirtley et al, Appl. Phys. Lett. 66, 1138 (1995). 
R. H. Koch et al, J. Low Temp. Phys. 51, 207 (1983). 
I. K. Kominis et al, Nature 422, 596 (2003). 
P. Treutlein et al, Phys. Rev. Lett. 92, 203005 (2004). 
M. Higbie et al, Phys. Rev. Lett. 95, 050401 (2005). 
B. Jo et al, cond-mat/0608585 (2006). 
L. Ho, Phys. Rev. Lett. 81, 742 (1998). 
Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 
(1998). 

R. H. Koch, D. J. Van Harlingen and J. Clarke, Appl. 
Phys. Lett. 38, 380 (1981). 

T. S. Lee, E. Dantsker and J. Clarke, Rev. Sci. Inst. 67, 
4208 (1996). 

Based on the geometry of the apparatus and the loca- 
tion of magnetic field sources, we estimate field variations 
across the condensate due to fourth order and higher 
terms to be <C 1 fT, negligible even at the high preci- 
sion of the spinor magnetometer. 

J. Dalibard and C. Cohen-Tannoudji, J. Opt. Soc. Am. 
B 6, 2023 (1989). 

H. Carr and E. Purcell, Phys. Rev. 94, 630 (1954). 



5 



[18] S. Wildormuth et al, Nature 435, 440 (2005). 
[19] A. Kuzmich et al, Europhys. Lett. 42, 481 (1998). 



